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Abstract 

We introduce an extension of the lattice model of melting of semiflexi- 
ble polymers originally proposed by Flory. Along with a bending penalty 
e, present in the original model and involving three sites of the lattice, 
we introduce an interaction energy Sp, corresponding to the presence of 
a pair of parallel bonds and an interaction energy Sh, associated with a 
hairpin turn. Both these new terms represent four-site interactions. The 
model is solved exactly on a Husimi cactus, which approximates a square 
lattice. We study the phase diagram of the system as a function of the 
energies. For a proper choice of the interaction energies, the model ex- 
hibits a first-order melting transition between a liquid and a crystalline 
phase at a temperature Tm- The continuation of the liquid phase be- 
low Tm gives rise to a supercooled liquid, which turns continuously into a 
new low-temperature phase, called metastable liquid, at Tmc < Tm- This 
liquid-liquid transition seems to have some features that are characteristic 
of the critical transition predicted by the mode-coupling theory. The exact 
calculation provides a thermodynamic justification for the entropy crisis 
(entropy becoming negative), generally known as the Kauzmann paradox, 
caused by the rapid drop of the entropy near the Kauzmann tempera- 
ture. It occurs not in the supercooled liquid, but in the metastable liquid 
phase since its Helmholtz free energy equals the absolute zero equilibrium 
free energy at a positive temperature. A continuous ideal glass transition 
occurs to avoid the crisis when the metastable liquid entropy, and not 
the excess entropy, goes to zero. The melting transition in the original 
Flory model, corresponding to the vanishing of the four-site interactions, 
appears as a tricritical point of the new model. 



1 INTRODUCTION 

Flory's model of crystallization of semiflexible linear polymers [1,2] is well known 
and describes a completely inactive crystal (CR) at low temperatures [3] within 
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the approximations developed independently by Flory [2], and by Huggins [4]. 
There is a first-order melting to a disordered equilibrium liquid (EL) phase, 
which has made it an attractive model to justify the Kauzmann paradox [5] and 
the ideal glass transition [6] in supercooled liquids (SCL) [5-31]. The ideal glass 
transition in fragile supercooled liquids [7,9,17] is believed to be a manifesta- 
tion of a rapid drop in the (configurational) entropy [32] near the Kauzmann 
temperature Tk: and has been a topic of many recent investigations [24-31]. 
The entropy in this work will always refer to the configurational entropy [32]. 
There are competing theories, both for and against an ideal glass transition, and 
the situation is far from clear [7]. Even the nature of the melting transition in 
the Flory model is in dispute [10-13,27], mainly because the Gujrati-Goldstein 
excitations [3,10,11] destroy the complete inactive nature of the crystal phase. 
The present work is motivated by this confused state of the field, and provides 
a convincing argument in favor of an ideal glass transition at a finite non-zero 
temperature. In order to substantiate our claims, we need to consider an exten- 
sion of the original Flory model of melting. We also clarify the nature of the 
melting transition in the Flory model [10,11,13,27]. Our conclusions are based 
on exact calculations. Some of the preliminary results have appeared earlier 
[28]. The present work provides the missing details in Ref. 28. 

According to the paradox, originally introduced by Kauzmann [5], the ex- 
trapolated entropy Ssct.{T) of the supercooled liquid becomes less than the 
entropy Scr{T) of the corresponding CR at temperatures below the glass tran- 
sition temperature Tq. It is a common belief [5-7,14,15] that the entropy of 
SCL cannot be less than the entropy of CR. However, it is worth noting [7,28,31] 
that there is no thermodynamic requirement for or against this. It also docs not 
violate the third law of thermodynamics. However, treating this possibility as a 
paradox, now conventionally known as the Kauzmann paradox or catastrophe, 
Kauzmann suggested that the system would either crystallize spontaneously [5] 
or undergo an ideal glass transition [5-10,17] to avoid the paradox. 

The existence of a glass transition caused by the above paradox has been 
originally justified [8] only in the Flory model of melting applied to linear poly- 
mers that are long. The approximate calculation [8] shows that the CR phase 
is completely inactive (zero specific heat, and zero entropy). The supercooled 
liquid avoids the Kauzmann catastrophe by undergoing a continuous transition 
called the ideal glass transition. This pivotal work enshrined the Kauzmann 
catastrophe as probably the most important mechanism behind the glass tran- 
sition. 

It should be stressed that the glass transition is ubiquitous and is also seen 
in small molecules. However, no model calculation exists that demonstrates the 
paradox for small molecules. Unfortunately, the approximations used by Gibbs 
and DiMarzio [8] have subsequently been rigorously proven to be incorrect, and 
unreliable [10-13], casting doubts on their primary conclusion of the existence of 
the Kauzmann paradox. Thus, there is currently no justification for the paradox 
as the root cause for the ideal glass transition, at least in long polymers. 

The current work is motivated by a desire to see if we can, nevertheless, 
justify a thermodynamic basis of the ideal glass transition in very long linear 



2 



polymers. To this end, wc perform an exact calculation. We should point 
out that recently we have discovered the existence of an ideal glass transition 
in a model of simple fluids [31(a)] and of dimers [31(b)]. However, this work 
deals only with long polymers. As discussed elsewhere [28, 31], the ideal glass 
transition in our view comes about not due to the originally proposed Kauzmann 
paradox caused by Sscl{T) <Scr{T), but because of the entropy crisis when 
the entropy of the state becomes negative. A negative entropy implies that 
there cannot be any realizable configuration of the system, which is impossible 
as there must be at least one configuration for the system to exist in Nature. 
Thus, in the following, wc interpret the Kauzmann paradox not in the original 
sense, but in the sense of the above entropy crisis. 

In the Flory model, a polymer chain is assumed to consist of n equal seg- 
ments, each with the same size as the solvent molecule. Each site of the lattice 
is occupied by either a polymer segment or a solvent molecule, and the excluded 
volume effects are taken into account by requiring a site to be occupied only 
once, either by a solvent molecule or a polymer segment. Wc can also think 
of the solvent as representing voids in the system. The polymer chain occupies 
a contiguous sequence of all the lattice sites connected by polymer bonds. For 
concretcness and ease of discussion, wc take the lattice to be a square lattice, 
which approximates a tetrahedral lattice on which the model is supposed to be 
defined. Both lattices have the same coordination number q=4. At every site, 
the polymer chain can assume cither a trans conformation (the conformation is 
related to the state of two consecutive bonds), when the consecutive bonds are 
coUinear, or one of the two possible gauche conformations, when the polymer 
chain bends. For a scmiflcxiblc polymer chain, every gauche conformation has 
an energy penalty e compared to a trans conformation. We set the energy for a 
trans conformation to be zero. The total energy of the system in a configuration 
on a lattice of N sites is 

£ =N,s, (1) 

where Ag is the number of gauche conformations present in the configuration of 
the system. This interaction involves three consecutive molecules of the chain 
and is the only one considered in the Flory model. 

No interaction between non-consecutive portions of the same polymer chain 
or between different polymer chains is taken into account in the Flory model 
since, according to Flory [2], the crystallization of polymers is not due to in- 
termolecular interactions but due to internal ordering/disordering and excluded 
volume interactions. Both the Flory [2] and the Huggins [4] approximations 
predict that the (configurational) entropy S{g) of the polymer chain for a given 
fraction g = Ng/N of gauche bonds goes to zero at a critical value go (where 
go is 0.45 in the Flory approximation [2, 10(a)] and 0.227 in the Huggins ap- 
proximation [4, 10(b)]). Correspondingly, the predicted entropy of the system 
is zero for any g < go and gives rise to the inactive phase for g < go- The 
result of the calculation is shown schematically in Fig. 1. The system is in a 
disordered liquid phase EL at temperatures higher than the melting tempera- 
ture Tm (curve BM). At Tm, the system undergoes a first order transition to a 
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Figure 1: Free energy vs. temperature for a semiflexible polymer: (a) Flory's 
calculation; (b) Gujrati-Goldstein upper bound. 



completely inactive ordered CR, characterized by a zero free energy and a zero 
density g (portion MO). There is a discontinuity in g at Tm- The results due to 
Flory and Huggins are qualitatively similar; the main difference is that Flory's 
estimate of Tm is about four times higher than that due to Huggins [10(b)]. 
However, the simulations [12.27] strongly support the presence of the Gujrati- 
Goldstein excitations that destroy the inactive crystal at low temperatures, but 
the nature of the melting transition remains uncertain, which makes the mathe- 
matical extrapolation MA representing the supercooled liquid [8] questionable. 
In particular, it is not clear if the extrapolation of the exact result would give a 
non-zero temperature where S{T) would go to zero but where g > 0. 

Rigorous lower bounds on S{g) per particle (and hence upper bounds for the 
free energy) as a function of g have been obtained [10,11]. Gujrati and Goldstein 
were able to prove that the entropy per segment of the chain in the case of a 
single polymer chain that occupies all the sites of the lattice (the Hamilton walk 
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limit) satisfies 

5(5)>(|)ln(l-3). (2) 

Hence, S is positive for any value of g >0, as it surely must be, in contrast 
with the results obtained by Flory. Bounds are also available for the case of 
finite- length polymers [11]. The above bound (2) implies that the equilibrium 
free energy of the system is never zero for T > 0, sec curve (b) in Fig. 1, and 
that the system is never completely ordered at any finite temperature. 

While the results due to Gujrati and Goldstein clearly show that the ap- 
proximations of Flory and Huggins do not give a satisfactory description of the 
system, they just provide an upper bound for the equilibrium free energy of the 
system; nothing is known about the correct equilibrium entropy. Therefore, it 
is still unknown what the actual behavior of the free energy is, whic;h is needed 
to obtain the continuation of the SCL liquid phase. The bounds do not say 
anything about the extrapolated (i.e. continued) SCL free energy or entropy. 
The knowledge of the reliable free energy form is fundamental in order to un- 
derstand if there is a phase transition of any kind in the system at any finite 
temperature. In particular, it is not clear if the model has a first-order melting 
transition. It should be recalled that one can usually continue the free energy 
only past a first-order transition, and not a second order transition due to the 
singularity in the latter case. If there is no melting transition with a latent 
heat, then there may be no SCL liquid below the melting transition. In this 
case, there would be no validity to the Gibbs-DiMarzio conjecture of a glass 
transition in the SCL liquid. Thus, an exact calculation is highly desirable. 

In recent years, the study of the glass transition has been stimulated by 
the development of the mode-coupling theory (MCT) [20-22] . This theory was 
developed in the first place for simple liquids but has been applied to polymers 
also [20]. The MCT studies the evolution of the density autocorrelation func- 
tion that can be measured in scattering experiments or calculated in computer 
simulations and is, therefore, of practical interest. The main result of this the- 
ory is the prediction of a critical temperature Tmc, above the glass transition 
temperature, corresponding to a crossover in the dynamics of the system. At 
Tmc, the correlation time of the system (the segmental relaxation time in the 
case of a polymeric glass) diverges with a power law just as one observes near a 
critical point: 

roc(T-TMc)"^. (3) 

Many neutron and light scattering experiments [20] have shown that the 
MCT is able to predict at least qualitatively the spectra of low molecular weight 
materials. Most of the systems for which MCT gives a good description of the 
dynamics (at least qualitatively) belong to the class of fragile glass formers. The 
theory has not been tested extensively with polymers that have large molecular 
weight but at the same time have been shown to be the most fragile systems 
yet identified [21-23]. 

Recent activities [24-26] liavc^ tried to export the progresses made in the field 
of spin glasses [33] to the field of real glasses. Even though the replica trick is 
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clearly unphysical [34,35], this approach has been extended to the case of real 
glasses. The replica approach has been applied to many Lennard- Jones glasses 
and the results have been interesting [24-26]. They provide some justification 
for the ideal glass transition. This may also be the case for polymers, which is 
the focus of this study. 

Despite the wide interest in the subject, there is still no comprehensive 
understanding of the nature of the vitrifying SCL and its relationship with CR, 
the mechanism responsible for the rapid entropy loss near Tq, and the nature 
of the ideal glass transition. It would also be interesting to see if there is a 
possible thermodynamic basis for the critical (and apparently a mode-coupling) 
transition in SCL's. 

In order to obtain a thermodynamic justification for all these phenomena, we 
consider in detail in this work a very simple limiting case. The solvent density 
will be taken to be identically zero. Thus, we consider an incompressible pure 
system. The effect of free volume is treated in separate publications [29-31]. We 
also consider the limiting ease of a single chain covering the entire lattice. Such 
a limit is conventionally known as the Hamilton walk limit [10,11]. The case of 
many chains of finite lengths is considered elsewhere [29-31]. To obtain a first 
order melting, we have to extend the Flory model of melting, as described below. 
We have substituted the original square lattice with a Husimi cactus (Fig. 2) 
on which the original problem is solved exactly. This is the only approximation 
we make. The results of this calculation for the case of a special interaction 
have been reported earlier [28] but details were not given. The present work 
also provides the missing details. 

It has been previously shown [36] that the exact calculations on recursive 
structures like the Husimi cactus are more reliable than conventional mean- 
field calculations. In this approach, the problem is solved exactly, taking into 
account all correlations present on the recursive lattice. In most cases, the real 
lattice is approximated by a tree structure. Because of the tree nature, the 
correlations are weak. We have chosen the Husimi cactus, obtained by joining 
two squares (Fig. 2) at each vertex, so that the coordination number 9 = 4. On 
a square lattice, there are also squares that share a bond. Such squares are not 
present in the cactus. Thus, the cactus should be thought of representing the 
checkerboard version of the square lattice, with the further provision that no 
closed loops of size larger than the elementary square are present. The square 
cactus is chosen to allow for the Gujrati-Goldstein excitations [10,11] that are 
important in disordering the ideal crystal at absolute zero. The results from the 
cactus calculation can be thought of as representing an approximate theory of 
the model on a square lattice. 

The layout of the paper is as follows. In the next section, we introduce 
the lattice model in terms of independent extensive quantities of interest. It 
is the most general model provided we restrict these quantities to represent 
pairs, triplets, and quadruplets of sites within each square. We also discuss 
the general physics of the model. As said above, we use a square lattice for 
simplicity to introduce the model, even though we eventually consider a Husimi 
cactus, on which the calculations are exactly. In Sect. HI, we explain the 
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Figure 2: Upper half of a Husimi cactus of generation m = 3. The dangling 
bonds outside it show its connection through surface sites (not shown in the 
figure) with the larger infinite cactus, as explained in the text. 



recursive sohition method on a Husimi cactus. We introduce 1-cyclc and 2- 
cycle solutions, representing the disordered and the ordered phase, respectively. 
The results are presented in Sect. IV along with the discussion, and the final 
section contains our conclusions. 

2 THE MODEL AND ITS PHYSICS 
2.1 Independent extensive quantities 

We consider a square lattice of N sites to focus our attention. We will neglect 
surface effects. There are Nb =2N lattice bonds, or distinct pairs of sites. Let 
us describe the state of a square by the number of polymer bonds j in it. The 
bonds in the following refer to the polymer bonds. Let A^'so, and Nsi denote 
the number of squares (S) with j = 0, and 1, respectively; see Fig. 3. 

For j = 2, we distinguish between the case of parallel bonds (p), with the 
number of such squares A^spj and gauche (g) bonds, with the number of such 
squares iVsg. The hairpin (h) turn corresponds to j = 3, with the number of 
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(c) » 



Figure 3: The possible states of a square in the lattice: (a) no bonds, (b) one 
bond, (c) two bonds and (d) three bonds. 



such squares A^sh- No square is allowed to have four bonds in it. Let Nt and 
A^g denote the number of trans and gauche bonds, respectively, and A^p and A^h 
the number of pairs of parallel bonds and hairpin turns, respectively, in a given 
configuration. We will also use them to represent their average values, as there 
will be no confusion. It is easily seen that the number of squares on a square 
lattice is A^s = AT. Let B denote the number of polymer bonds, and A^mm the 
number of unbonded monomer- monomer contacts. The following topological 
identities are easily seen to hold: 

A^s = A^so + A^si + A^sp + A^sg + A^sh, (4) 
2Afmm = 4Afso + SNsi + 2Nsp + 2Nsg + ATgh, (5) 
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2B = Nsi + 2iVsp + 2Nss + SA^sh, 



(6) 



N = Nt + N^, 2N = B + N^^, (7) 

A^p = iVsp + A^Sh, A^g = TVsg + 27Vsh. (8) 

As said earlier, the cactus represents the checkerboard version of the square 
lattice, so that the number of squares Ns on both lattices with the same number 
of sites A'' are not the same; see Sect. Ill also. For a square lattice, -A^s= A/'; 
for the cactus, Ns = N/2. However, the number of lattice bonds Nb on both 
lattices are the same. Because of this, Eq. (6) must be modified for the cactus. 
Since each bond belongs to only one square in the cactus, we have: 

B = Nsi + 2Nsp + 27Vsg + SiVgh. (6a) 
All other identities remain valid on the cactus. 

2.2 General Model 

Among the eleven extensive quantities, N, B, TVg, A'^t, A^mm, A^so, A^si, A^Sp) 
A^Sg) A'^st and A^sh, there are six independent geometrical relations; the second 
one in (7) is not independent. In addition, for the Hamilton walk problem, 
we have B = N. Thus, there are only four independent extensive quantities, 
which we take to be N, Ng, Np and A^h. One of these, the lattice size A'', 
will be used to define the partition function. The remaining three independent 
quantities A'p, A'g and A^ii will bo then used to define the configuration uniquely. 
Corresponding to each of the quantities A^g, A^p and A^hi there is an independent 
activity w, Wp and w^, respectively, that will determine the partition function 
for the Hamilton walk problem as 

^iv = 5^.«^«<''<N (9a) 

where the sum is over distinct configurations obtained by all possible values of 
A^g, Np and A^h consistent with a fixed lattice size N. The activities w, Wp 
and uih arc determined by the three-site bending penalty e introduced by Flory 
in his model, an energy of interaction Sp associated with each parallel pair of 
neighboring bonds, and an energy Sh for each hairpin turn within each square 
as follows: 

w = exp (-/?£) , Wp = exp (-/Jffp) , Wh = exp{-l3sh). 

Here, /3 is the inverse temperature T in the units of the Boltzmann constant. 
The original Flory model is obtained when the last two interactions are absent. 
It should be stressed that e\i is the excess energy associated to the configuration, 
once the energy of the two bends and the pair of parallel bonds have been 
subtracted out. Both Sp and £h are associated with four-site interactions, since 
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it is necessary to determine the state of four adjacent sites to determine if a pair 
of parallel bonds or a hairpin turn is present. 

The model is easily generalized to include free volume by introducing voids, 
each of which occupies a site of the lattice. The number of voids A^o is controlled 
by the void activity 77. The number P of polymers is controlled by another 
activity given by H^. The interaction between nearest neighbor pairs Nc of 
voids and the monomers of the polymers determines the Boltzmann weight Wq. 
The partition function of the extended model is given by 

Zm = Y^ J?'^w;^«;^«<-<N (9b) 

where the sum is over distinct configurations consistent with the fixed lattice 
of N sites. Because the activity H only determines the average number of 
linear polymers, but not their individual sizes, the model in Eq. (9b) describes 
polydisperse polymers [37], each of which must contain at least one bond. 

We now turn to our simplified model of the Hamilton walk (P = 1, and 
A^o = 0). In this model, the energy of interaction in a given configuration is 
given by 

£ = eN^ + SpNp + EhA^h = e {Ng + aNp + feiVh) , (10) 

where a = Sp/e, b = Sh/e. The parameters can in principle assume positive and 
negative values. However, we will restrict ourselves to e > in this work. The 
limit e = Sp = Sh = corresponds to a completely flexible polymer problem, 
which is of no interest to us here, as it corresponds to the infinite temperature 
limit of our model. The limit, however, is of considerable interest in the study of 
protein folding and has been investigated by several workers [38]. In addition, 
we will focus mainly on the case < a < 1. A positive a guarantees that 
parallel bond energy opposes the creation of configurations in which pairs of 
parallel bonds arc present and a < 1 makes the penalty for a pair of parallel 
bonds less than that for a gauche conformation. This guarantees the presence 
of a crystalline phase at low temperatures, as shown below. 

2.3 Ground State at T=0 

The physics of the model at absolute zero can be easily understood on general 
grounds. We are interested in the thermodynamic limit N 00. We first 
consider 6 = 0. For a < 1, the ground state at T = has Ng = 0, Np = N, 
= 0, as shown in Fig. 4(a). (The labels R and L are related to the state 
of the sites as introduced in the next section.) Thus, £ = SpN. This is what 
we will call the perfect crystal at absolute zero. For b ^ 0, the state in Fig. 
4(a) remains the ground state as long as 2 + & > 0. This condition ensures that 
hairpin turns are not present. For a > 1, and 6 = 0, the ground state at T = 
has Ng = N, Np = 0, as shown in Fig. 4(b). Thus, £ = eN. This remains the 
ground state provided a + 6 > 0, which ensures that the hairpin turns are not 
present. Since our interest is to have a crystal state as the equilibrium state at 
low temperatures, we would only consider the earlier case a < 1, with 2 + 6 > 0. 
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Figure 4: Possible configurations of the polymer chain at T = 0: (a) crystalline 
phase with iVg = and A^p = iV; (b) step-like configuration with iVp = and 
A^g = N. See text for the explanation of the symbols. 



It should be recognized that the model considered here is defined on a lat- 
tice. Thus, the ground state also possesses the symmetry of the lattice. This 
symmetry is imposed by the lattice symmetry and is independent of the model. 
Thus, this induced symmetry should not be confused with the point group sym- 
metry of a real crystal, which is brought about by the interactions in the system. 
The symmetry in our model is due to the orientational order between pairs of 
parallel bonds. It is because of the orientational order defining the crystal in 
our model that we can obtain a continuous transition between the crystal and 
the equilibrium liquid. 

The lack of a point group symmetry of a real crystal in our model should not 
be a taken as a serious limitation of the model, since our main goal is to study 
the possibility of a glass transition in a supercooled liquid. The determination 
of the supercooled liquid branch requires the continuation past a first-order 
melting transition. Thus, the exact nature of the symmetry of the CR phase is 
not as important as the existence of a discontinuous melting. 
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3 RECURSIVE SOLUTION 



The Husimi cactus approximates the square lattice, as said above. Both have 
the coordination number q = A, and tlie elemental squares as the smallest loop. 
However, the most important reason for choosing the square cactus is that it 
allows for hairpin turns that give rise to the Gujrati-Goldstein excitation in 
the Flory model of melting. A Bethe lattice would be inappropriate for this 
reason. As said above, the number of squares on the cactus is half of that on 
a square lattice with the same number of sites A'^. This can also be easily seen 
as follows by assuming homogeneity of the lattice. First, consider the square 
lattice. Four squares meet at each site; however, each square will be counted 
four times, due to its four corners, assuming homogeneity. Thus, A^s = N. On 
a cactus, only two squares meet, but each one is counted four times as before. 
Hence, A's — N/2. Despite this, A^b = 2A^ on both lattices, only half of which 
are going to be taken up by the Hamilton walk on both lattices. 

A site is shared by four bonds and two squares S, and S' that are across 
from each other on the cactus. On the other hand, there are two different pairs 
of such squares on a square lattice. In a formal sense, we can imagine that each 
end of a bond contributes | of a site, and each corner of a square contributes 5 
of a site (on a square lattice each corner contributes ^ of a site.). This formal 
picture will be useful in determining the nature of a homogeneous cactus. To 
make the cactus homogeneous, we must consider it to be part of a larger cactus. 
This is shown in Fig. 2, where we show a cactus of generation m = 3 with 
dangling bonds (each one ending with a surface site, not shown in the figure) 
outside it to show its connection with the larger infinite cactus. The latter has 
no boundary. A similar homogeneity hypothesis associated with a Bethe lattice 
has been discussed in Ref. 37(b) to which we refer the reader for further details. 
On a Bethe lattice, each dangling bond was treated as a half-bond to ensure 
that A^B = qN/2, where q is the coordination number of the Bethe lattice. For 
the case of the cactus, we treat each pair of dangling bonds in Fig. 2 as a 
half-square, and each surface site as a half-site to calculate the number of sites 
Nm, and the number of squares Sm for a cactus of m generations. A trivial 
calculation shows that 

Ar„ = 4 X 3™ - 2, 5„ = 2 X 3™ - 1, (11) 

so that Sm = Nm/2, as m — > 00. Since each square contributes 4 bonds, it 

is also evident that A^b = 2A'' in the limit of an infinite cactus. A detailed 
calculation of the quantities introduced above is given in the Appendix. 

Because of the above-mentioned homogeneity, a site is arbitrarily designated 
as the origin of the cactus. Each square has one site, called the base site, closer 
to the origin. The base site is given an index m > , the two sites next to 
the base site within the square, called the intermediate sites, the index (m-|- 1), 
and the remaining fourth site, called the peak site, the index (m + 2). We 
will call this square an mth level square; it has its base at the mth level and 
its peak at the (m + 2)th level; see Fig. 5. The two lower bonds in the mth 
square connected to the mth site are called the lower bonds and the two upper 
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Figure 5: The four possible states of the polymer chain at any site at the mth 
level of the lattice. 



bonds connected to the peak site are called the upper bonds. The origin of the 
lattice is labeled as the m = level and the level index m increases as we move 
outwards from the origin. We can imagine cutting the Husimi tree at an mth 
site into two parts, one of which does not contain the origin if m > 0. We call 
this the mth branch of the lattice and denote it by Cm- At the origin, wc get 
two identical branches each containing the origin. We will call each of those the 
(m — 0)th branch. 

We only consider parallel bonds and hairpin turns that are inside the squares, 
since the cactus represents the checkerboard version of the square lattice. Thus, 
each square can contribute only once to either A^p, or A^h- This means that 
Np = Ns = N/2 in the perfect crystal at absolute zero. Similarly, it can be 
easily seen that the maximum possible value of A^h is 2iVs/3 = iV/3. 
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3.1 Recursion Relations 

Wc consider a linear polymer that covers all the sites of the Husimi cactus or 
the square. Its configuration determines the state of the bonds in each square. 
Consider a pair of two squares E, and S'that are across from each other. We 
distinguish E by putting a filled dot (•) just above the common site. We now 
face towards E' from within E through this common site. The common site has 
been taken as the base site in Fig. 5, but the following description is valid at 
any site. The common site can assume 4 possible different states depending on 
the state of the four bonds connected to it. Two of the bonds are in E, and the 
remaining two are in E'. 

1. In the I state, both E- bonds are occupied by the polymer chain. Since the 
polymer is linear, the two E'-bonds must be unoccupied by the polymer. 

2. In the O state, both E-bonds are unoccupied but both E'-bonds are oc- 
cupied by the polymer chain. 

3. In the L state, only one of the E-bonds is occupied and the polymer 
occupies the left bond in E' (we always think about left and right as we 
face towards E'). 

4. In the R state, only one of the E-bonds is occupied and the polymer 
occupies the right bond in E'. 

For the common site at m = 0, the square E' in the above classification is 
the square on the other side of the origin. 

It is now easy to understand the labeling of the two configurations shown in 
Fig. 4. 

We are interested in the contribution of the portion of the mth branch C,„ 
of the lattice to the total partition function of the system. This contribution is 
called the partial partition function (PPF) of the branch. It is easy to see that 
the PPF depends on the state a of the mth lev{4 site. We denote this PPF by 
Zm{ce). We now wish to express Zm{a) recursively in terms of the PPF's of the 
two intermediate sites and the peak site of the roth square. Following Gujrati 
[36] , the recursion relations can always be written in the form: 

Z„(a) = Tr[W(a,{/3})Z^+i(/3i) mzm+2m], (12) 

where {/?} is the set of states (3^, the latter states represent the possible states of 
the other three sites of the square, and W(a, {/?}) is the local Boltzmann weight 
of the square due to conformation of the polymer chain inside the square. 

Let us consider in detail the case in which the base site at the roth level is in 
the I state. The three possible configurations that the polymer can assume in 
the roth level square are shown in Fig. 6. In this figure, L, R, I and O represent 
the possible state of the (ro -|- l)th and (to -I- 2)th level sites and w represents 
the weight of a bend. In order to carefully account for statistical weights, a 
Boltzmann weight equal to w is considered only if the bend happens: 
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(1) 



Figure 6: Possible configurations of the polymer chain when the mth level site 
is in the I state. 



1. At the mth level and at least one polymer bond at the level is inside the 
square. 

2. At the (m + l)th or (m + 2)th level, and both polymer bonds at the level 

are inside the polymer. 

A weight Wp = w"" is considered for any configuration in which two bonds are 
parallel to each other within the same square. We can, furthermore, distinguish 
configurations in which two bonds are parallel to each other from configurations, 
like those shown in Figs. 6(b)-(c), where three consecutive bonds form a hairpin 
configuration. Whenever this configuration is present, an additional weight 
w\i = is introduced. 

It is not important to know akmg which of the two lower bonds in the 
(m + l)th or (m + 2)th square docs the polymer chain enter into the mth 
square. In fact, even if the polymer undergoes a bend while moving from the 
higher level square to the mth level square, the corresponding weight is already 
taken into account into the partial partition function of the higher level site. 

It is important to consider always the state of a site moving through the 
lattice towards the origin. In the configuration in Fig. 6(a), the intermediate 
site on the left is in the R state since the polymer undergoes a right turn after 
entering the square. The intermediate site on the right is in the L state since 
the polymer undergoes a left turn after entering the square. Finally, the peak 
site is in the I state because the polymer is coming from the (m + l)th level 
square but not entering the mth level square. The polymer undergoes one bend 
(at the base site) so that there is a weight w to be taken into account. Thus 
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the contribution to the partition function coming from this configuration is: 

wZm+l{R)Zm+l{L)Zm+2i'i)- (13) 

In the configuration in Fig. 6(b), the intermediate site on the left is in 
the O state since both the lower bonds in the corresponding (m + l)th square 
are unoccupied. The intermediate site on the right is in the L state since the 
polymer undergoes a left turn after entering the square. Finally, the peak site is 
in the R state since the polymer undergoes a right turn after entering the square. 
The polymer undergoes two bends (one at the base site and the other one at 
the left intermediate site) so that there is a weight to be taken into account. 
There is a pair of parallel bonds in the square and a hairpin turn occurs so that 
a weight WhWp has also to be taken into account. Thus, the contribution to the 
partition function coming from this configuration is: 

wWwpZm+liO)Zm+imZm+2{R). (14) 

In the configuration in Fig. 6(c), the intermediate site on the right is in 

the O state since both the lower bonds in this particular (m + l)th square are 
unoccupied. The intermediate site on the left is in the R state since the polymer 
undergoes a right turn after entering the square. Finally, the peak site is in the 
L state since the polymer undergoes a left turn after entering the square. The 
polymer undergoes two bends (one at the base site and the other one at the 
right intermediate site) so that there is a weight to be taken into account. 
We also have a pair of parallel bonds and a hairpin turn to take into account 
in this case. Thus, the contribution to the partition function coming from this 
configuration is: 

w'^WhWpZm+l{0)Zm+l{R)Zm+2m- (15) 

The recursion relation for Zm(i), the partition function of the mth branch 
of the Husimi tree given that the mth level site is in the I state, is therefore 
given by: 

Zm{l) = w'^Wi,WpZ,n+l{0)[Z„i+l(L)Zm+2(R) + (R)-^m+2 (L)] 

+ w;Z„+i(R)Z„+i(L)Z„+2(I)- (16) 

Considering the case in which the mth level site is in the O state, the partial 
partition function Zm{0) for the O state can be written as: 

■^m(O) = Z^_^;^(I)Zj„+2(I) + Zjn+l{l) [Zm+l(J^)Zm+2(R) + ■^m+l(R)^m+2(L)] 
+ wZm+l{R)Zm+l{L)Zm+2{0). (17) 

When the mth level site is in the L state the partial partition function can 
be written as: 

■^m(L) = [Zm+l{R) +wZm+l (L)]{Z„+i(I)Z„+2(I) + W^WpWhZm+l{0)Zm+2{0)} 

+ Wp [wZl+,{L)Z„,+2iR) + ^^+i(R)^™+2(L)] 

+ [^„+2(R) + wZm+2M wZm+l{I)Zm+l{0). (18) 
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The relation for the R state is obtained from Zm{h) by the interchange 

L< — >R: 

■^m(R) = [2m+l(L) + wZ„i+i(R)]{2'„i+i(l)2'„i+2(I) + W^W-pWh.Zm+l{0)Zm+2{0)) 

+ [Zm+2{L) + wZ^+2iR)]wZ„,+iiI)Z^+i{0). (19) 

It is possible to write analogous relations for Zm+i{ce) by properly substitut- 
ing m— >m+l,m+l— >TO + 2 and to + 2 ^ m + 3. We introduce the following 
ratios between partial partition functions at even and odd levels of the lattice: 

Xm{l) = Zm{T)/ {Zm{L) + ^m(R)) , 
Xm{0) = Z^(Qi)l (Z„(L) + Z„(R)) , 

a;„(L) = Z,„(L)/ (Z™(L) + Z„(R)) , 

a;„(R) = l-a;„(L). (20) 

As one moves from a level that is infinitely far away from the origin towards 
the origin itself, the recursion relations (16)-(19) will approach fix-point (FP) 
solutions, Xm{oi) x*{a), Xm+i{a) x**{a), etc., where a = I, O, L or 
R. These fix-point solutions of the; recursion relations describe the behavior in 
the interior of the Husimi tree. Once the fixed point is reached, the value of 
x*{a) and x**{a) becomes independent of to. On a Husimi cactus, a site can 
be classified as a simultaneous peak and a base site, or a simiiltaneous peak 
and a middle site, depending on the pair of squares which share the site. Thus, 
it is expected that the most general FP solutions will correspond to a 2-cycle 
solution in which Xm{a) and Xm+2{<^) tend to the same limit. In this case, we 
obtain a sublattice structure in which sites with even levels are different from 
sites with odd levels. We can write in this case: 

Xm(i) = Xm+2{1) = «a, ^^^(O) = Xm+2iO) = Oa, ^^(L) = Xm+2(L) = /a, 
Xm{R.) = Xm+2(R) = 1 - ^^(L) = 1 - l^, 
Xm+li^) = Xm+si^) = ib, a;m-|-l(0) = Xm+3{0) = Ob, Xm+l^ = a;m+3(L) = Zb, 
Xm+l{R) = Xm+3{^) = 1 - Xm+l{l-) = 1 - Ih- (21) 

The indices a and b refer to even and odd levels, respectively. Using (21), it 
is easy to prove that the system of equations (16) to (19) can be written in the 
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following form: 

«aQLR = W{1- lb) l^i^ + W^WpWhOb [/b(l - ^a) + ^a(l " ^b)] , (22) 

OaQLR = «b^a + ih [hO- ' Q + ^a(l - ^b)] + w{l - h)ho^, (23) 
^aQLR = (1 - /b + w/b)[«b«a + W.'^Wp WhOfaOa] 

+ 10Obib(l - ^a + wis) + Wp(l - Ihfh + WW^ll{l - la), (24) 

^b^LR = W (1 - D IJh + W^WpWhOa [^a(l " ^b) + ^b(l " ^a)] , (25) 

Ob^LR = «a«b + ia [^a(l - ^b) + ?b(l - ^a)] + w{l - k)kOb, (26) 
^bQLR = (1 - ^a + Wk)[ijh + W^U'pWhOaOb] 

+ WOaia(l - ^b + Wlb) + Wp{l - Za)^^b + WWpll{l - lb), (27) 

whore is obtained from Qlr by exchanging a and b subscripts and Qlr 

can be written as: 

QlR= (1 + w) {ibia + W^U;pWhObOa + wibOb + Wp[Zb(l-^a)^+^a(l-^b)^]}- (28) 

3.2 2-cycle Free Energy 

In order to determine which phase is the stable one at some temperature, we 
have to find the free energy of all the possible phases of the system as a function 
of w. We follow the treatment by Gujrati [36], and provide its trivial extension 
to the 2-cycle FP solution shown above. The free energy per site at the origin of 
the lattice can be easily calculated from the expressions for the total partition 
function Z at the (m = 0)th, (m = l)th and (m = 2)th levels. 

The total partition function of the system Zq can be written by considering 
the two (m = 0)tli branches Co meeting at the origin at the (m = 0)th level. 
For this, we need to consider all the possible configurations in the two branches. 
This is done by considering all the configurations that the polymer chain can 
assume in the two squares that meet at the origin of the tree. All the possible 
configurations of the (m = 0)th level site (in this case we are not interested in 
the state of the (m = l)th level sites) are shown in Fig. 7. Each of the first two 
configurations contributes 

Zo{l)Zo{0) 

to the total partition function. The third and fourth configurations both con- 
tribute 

(l/«;)Zo,g(L)Zo,g(R), 

where the factor (1/w) is needed in order not to take into account the Boltzmann 
weight at the origin twice and the subscripts "g" and "t" refer to the gauche 
and trans part of the partition function for L and R states. In fact, it is possible 
to separate these two contributions to the partition function at any level. The 
"gauche" portion is the one corresponding to configurations such that there is a 
bending at the mth level site, while the "trans" portion is the one corresponding 
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Figure 7: Possible configurations of the polymer chain at the origin of the tree. 



to configurations in which the two bonds coming out of the mth level site that 
we are considering are straight. It is easily seen that: 

■^m,t(L) = Zm+l(R){Z„i+i{l)Zrn+2(i) + (0)Z„+2 (O) 

+ WpZ„+i(R)Z„+2(L)} +m;Z„+2(R)Z™+i(I)Z„+i(0), (29) 

and 

Zm,s(L) = W^m+l(L){Z„i+i(I)Zm+2(I) + w'^ WpWhZ„i+l{0) Zm+2{0) 

+ WpZ^+i{L)Z^+2{R)} + w^Z^+2{L)Z^+i{l)Z^+i{0). (30) 

Finally, the fifth and sixth configurations contribute Zq^{L) and .^q^(R), 
respectively. It is then possible to write: 

Zo = 2Zoil)Zo{0) + (2/u;)Zo,g(L)Zo,g(R) + Zo%(L) + Zl,{R). (31) 

It is clear that Zq is the total partition function of the system obtained by 
joining two branches Co together at the origin. Now, let us imagine taking away 
from the lattice the two squares that meet at the origin. This leaves behind four 
different branches Ci and two branches C2. We connect the two C2 branches to 
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form a smaller cactus whose partition function is denoted by Z^. Similarly, 
we join two of the C\ branches to form an intermediate cactus whose partition 
function is denoted by Z\ . We can form two such intermediate cacti out of the 
four Ci branches. Each partition function Z\ or Z^ can be written in a form 
that is identical to that of Equation (31): 

Z, = 2Zi{I)ZiiO) + (2/«;)Zi,g(L)Z,,g(R) + Zl{L) + Zl{R), (32) 

whore i = 1, 2. 

The difference between the free energy of the complete cactus and that of 
the three reduced cacti is just the free energy corresponding to a pair of squares 
so that, following Gujrati [36], we can write the adimensional free energy per 
square (without the conventional minus sign) as: 



o,q = -ln[Zo/{Z^Z2}]. (33) 



It is possible to write: 

Zq = BlQ2{ia., Oa, ^a, «b, Ob, ^b), (34) 

Zi = BlQ2{ia.,0a.,la.,ih,0h,lh), (35) 
Z2 = -B|(52(ia, Oa, la., ib, Ob, lb), (36) 

where we have introduced 

B,n = ZmiX) + Zm(J^, (37) 

and where Q2 is the following polynomial of ia, Oa, ^a, ^b, Ob and Zb- 

Q2 = 2iaOa + (2/«;)Za,g(l - /a)g + ^a.t + (1 " Qtl (38) 

/a,t and /a,g correspond to the trans and gauche portions of respectively, and 
Q2 is obtained from Q2 by interchanging a and b subscripts. 

It is also easily seen that: 

Bo = BlB2QhK{ia, Oa, ^a, «b, Ob, Ih)- (39) 

so that the free energy per square can be written as: 



In ( ^ ) . (40a) 



1^2 

The free energy per site Wgite is proportional to Wgq : 



sq 

Wsite = a;sq/2, (40b) 
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since there are two sites per square. 

The usual Helmholtz free energy per square can be obtained from uj through: 

F = -Tuj. (41) 

If it happens that the even and odd sites are not different, we obtain a 1-cycle 
FP-solution. Below, we will consider the two solutions separately. 

3.3 1-cycle solution 

In the 1-cyclc scheme, we have ^^(a) = Xm+i{oL) as they converge to the same 
fix point. Thus, we have: ia = *b = Oa = Ob = o and /a = = In this case 
the system of equations reduces to: 

iQLR = w{l-l)li + 2v?WpWY,ol{l - I), (42) 

oQlr = i^ + -l) + w{l - l)lo, (43) 

IQ-LK = [(1 - + wl] {i^ + w^WpwW + Wp«(l - + wio] , (44) 

with: 

<3lr = (1 + w) [i^ + vo^WpwW + Wp/(1 - + wio) . (45) 

From (44), it is easy to show that we must have Z = 5 for every solution obtained 
in this scheme. 

One solution that exists for every value of w (and, hence, of T) is the fol- 
lowing: 

/ = ^,o = 0,i = 0. (46) 

This represents a liquid-like phase without any I, or O states. We label this phase 
metastable liquid (ML) because, as we will show below, it never represents the 
equilibrium phase of the system even though it exists at all temperatures. For 
large enough w, there is another solution of the system of equations with nonzero 
values for o and i, and it has to be found numerically. We label this second 
liquid-like phase equilibrium liquid (EL) since, as our free energy calculations 
will show, it represents the equilibrium phase of the system at high temperature. 
The temperature at which EL appears is a function of Wp and w^. We will 
call the temperature Tmc, at which the equilibrium liquid appears, the mode 
coupling temperature for reasons that will become clear below. 

Table 1 shows how the value of Tmc changes as a function of a (both positive 
and negative values of a are considered, see below) as we keep b equal to zero. 
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a 


Tmc 


-1 


2.9586 


-0.8 


2.(il83 


-0.5 


2.1876 


-0.2 


1.7359 





1.4427 


0.2 


1.1600 


0.5 


0.7653 


0.8 


0.4156 


1 






Table 1: Values of Tmc as a function of a (with 6=0). 
For the ML, we have o = i = 0, Z = | so that: 

1 



and 



Q2 



2(1 + w)' 

Wp{l + w) 



so that the ML free energy per square assumes the simple form: 



WML = In 



Wp {l + wf /2 



In 



il + wy/2 +\n{wp), 



(47) 
(48) 

(49a) 



while, for the EL, we have to substitute the numerical solutions i{w), o(w), and 
I = I obtained from Eqs. (42)- (45) in the expression for the free energy. The 
ML entropy per square is given by 



5ml = In 



(1 



72 +2w/[T{l + w). 



(49b) 



The corresponding energy is given by 



-Bml = 2w/{l + w)+a. (49c) 

It is easily seen that the ML specific heat is always non-negative. It is very 
important to observe that the free energy of the ML phase does not depend on 
the value of the parameter a (except for the additive factor In(wp)) while the 
free energy of EL strongly does. At absolute zero, the ML entropy and energy 
go to ln(l/2) , and a, respectively. We find the modified free energy 

F = F-ep (50) 

more convenient to use than F itself since, at T = 0, the ground state is the 
one in which all the bonds are parallel to each other and the free energy of the 

system is equal to £p so that the crystalline ground state has always F = 0, 
regardless of the value of Sp. The free energy curves for the EL and ML phases 
are shown in Fig. 8. We immediately observe that ML at very low temperatures 
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Figure 8: a — 0.5, 6 = 0. Free energy in the 2-cycle FP scheme for the ML 
(continuous Une), EL/SCL (dashed hne) and CR (dash-dot hne). We also show 
7cre(4), Tm{o) and Tmc(»)- Here, as well as in Figures 10 to 12, the stable 
phases are represented by thick lines while the metastable phases are represented 
by thin lines. 



T < 0.48 (dash-double dot line in Fig. 8 originating at the origin) has negative 

entropy, since its free energy F is increasing with the temperature. A negative 
entropy is not possible for states that can exist in Nature. 

3.4 2-cycle solution 

The phase diagram obtained in the 1-cycle solution scheme cannot be complete 
because, at low temperatures, ML cannot be the stable phase. At T = 0, CR 
contains an alternating ordered sequence of L and R states in addition to having 
/ — i and no O and I states; see Fig. 4. This is a 2-cycle pattern in L and R 
that is completely missed by the 1-cycle calculation performed above. For ML, 
I = i also, but L and R are statistically distributed. One of these distributions 
must be the crystal state at T = 0; indeed, Fml(T = 0) = Fgr{T = 0). 
Despite this, ML immediately above T = can not represent CR, as it has 
negative entropy. To obtain the alternating sequence in CR at T > 0, the above 
1-cycle FP scheme is not sufficient to completely describe the physics of the 
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system. Wc also observe that for T > 0, there must be local Gujrati-Goldstcin 
excitations [10,11] creating imperfections by local L< — >R interchanges in the 
ordered [..LRLRLR..] sequence. The excitations change a local string LRL 
into LLL, or RLR into RRR within a square and require 4 bends only. Other 
excitations, which require (L or R)< — >(I or O) on the cactus, can not be done 
locally and require infinite amount of energy, and need not be considered. This 
means that the local density / or r will no longer be i. However, if Z > i at 
some site, then r > ^ at the next site, followed by I > ^ on the next site and so 
on. 

There are three solutions for the complete system of Eqs. (22)- (27) for any 
given value of the weights w, Wp and Wh- 

i. A metastable liquid ML (already found in the 1-cycle FP scheme) with: 

k = h = rg, = rb = ^ and ia = ih = Oa = Ob = 0. 

As seen above, this phase represents a liquid phase in which no O and I 
states are present. The R and L states are randomly distributed in the lattice 
with the only constraint of having the same number of L and R states at both 
odd and even layers. This solution exists for any temperature and its free energy 
has a maximum at T = Tk — 0.48. 

ii. An equilibrium liquid EL characterized by the presence of all possible 
states I, O, L and R at both odd and even levels, so that 

k = lh = ra. = rb = l and i^,, i^, o^, 7^ 0. 

In the 2-cycle solution, i and o on the two sublattices are different, which 
makes this solution different from the 1-cycle EL solution, in which there is no 
sublattice structure. Despite this, EL phases in both schemes have identical free 
energy and various densities. Thus, we no longer make any distinction between 
the two solutions and identify both of them as the same EL phase. As seen in 
the previous subsection, the free energy of this phase depends on the value of 
the parameters a and b. This solution exists only for temperatures larger than 
T = TMc(a, 6). 

iii. A crystal phase CR with double degeneracy that is the ground state 
and that exists for temperatures lower than T = Tcre = l/lii(2)- The state 
is perfectly ordered at zero temperature and disorders as the temperature is 
raised. Fig. 9 shows how the values of la and If, change with temperature for 
the CR phase: the two degenerate solutions correspond to a different labeling 
of the lattice sites where the odd and even levels are just exchanged with each 
other. 

The solutions of the system of Eqs. (22)-(27) corresponding to CR and 
to ML do not depend on the strength of the three- and four- site interaction 
and, therefore, the free energy curves corresponding to these two phases do not 
change when the parameters a and b change. In contrast, the free energy of the 
EL phase depends on the value of a and b. 

For a < 1, we have the ground state at T = in which all the bonds are 
parallel and there are no bends in the systems. If a > 1, the four-site interaction 
becomes more important than the three-site interaction and the ground state at 
T = is a step like walk in which there are no parallel bonds and all the sites 
of the lattice are in a gauche conformation (so that = Ns, as said earlier). 
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Figure 9: Dependence of l^, and on the temperature for the two phases ob- 
tained at low temperature in the 2-cycle fixed point scheme. 

The two possible ground states are shown in Fig. 4. 

4 RESULTS AND DISCUSSION 

4.1 Thermodynamic Functions 
4.1.1 b = 

The complete free energy diagram for a = 0.5 is given in Fig. 8. The equilib- 
rium phases are represented by the disordered EL at high temperatures and the 
ordered CR at low temperatures, with a first order transition at a temperature 
Tm between the two phases, with a discontinuity in the first derivative of the 
free energy with temperature. This remains true as long as a > 0. The situation 
with a < is different and is discussed later. 

The existence of a discontinuous melting temperature for o > makes it 
possible to have a supercooled liquid phase through continuation. For T > Tm, 
the EL phase is the stable one. If the liquid phase is cooled in such a way that 
it is not allowed to undergo the melting transition at Tm, then it is possible 
to have (for T < Tm) a supercooled liquid (SCL). The free energy of SCL 
is obtained by continuing the free energy of the EL phase. This free energy 
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meets critically (i. e. with continuous slopes) with the ML free energy at a 
temperature that, as before, we call Tmc- For a = 0.5 we find that Tmc > ^k, 
where Tk is the temperature where the ML free energy has its maximum. The 
critical transition between ML and SCL is a liquid-liquid transition between 
two liquid phases. As a increases, Tmc moves towards Tk- We have observed 
that for a > 0.8, Tmc < (results not shown). In particular, the EL free 
energy curve itself has a maximum in this case before it merges with the ML 
curve and, consequently, has an unphysical portion corresponding to the entropy 
crisis below its maximum. 

Figs. 10 and 11 show the entropy and the specific heat vs. temperature, 
respectively, corresponding to the free energy results shown in Fig. 8. As 
explained before in the case of the free energy-temperature graphs, the curves 
corresponding to CR and to ML do not depend on the choice of a and b. Table 
2 shows how the value of Tm changes as a function of a (in the case 6 = 0). 



a 


Tm 





1.443 


0.2 


L351 


0.5 


L198 


0.8 


1.009 


1 


0.878 



Table 2: Values of Tm as a function of a (with 6 = 0). 

Only positive values of a are considered since, as it will be shown below, 
when the parameter a is negative, there is no first-order melting in the system, 
provided that 6 = 0. 

We can calculate the density of gauche bonds g and the density of pairs of 
parallel bonds p as a function of T and a. We can write: 

g = du)site/d{\n{w))\yj^,i^^, (51) 

(52) 

and calculate the two densities from these derivatives. Note that we have defined 
the gauche bond density g per site, while the density of parallel bond pairs p is 
defined per square, since each square contributes one such pair in the ideal CR 
at absolute zero. Fig. 12 shows the gauche bond density and the parallel bond 
density in the case of a = 0.5, 6=0. 

4.1.2 6^0 

The effect of changing the value of the parameter 6 is shown in Fig. 13 for 
a = 0.5. As we can see, a change in 6 does not have any effect on the free energy 
of the CR and ML phases but it does affect the EL phase. Apparently, the effect 
of 6 is smaller than the effcict of n, since the value of the melting temperature 
does not change significantly as we change 6. It is worth noting that the presence 
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Figure 10: a = 0.5, 6 = 0. (a) Entropy in the 2-cycle FP scheme for the ML 
(continuous Une), EL/SCL (dashed hne) and CR (dash-dot hne). We also show 
Tcre(4), Tm{o) and Tmc{*)', (b) magnification of the area contained in the box 
in (a). 



of the hairpin term alone, even in the absence of the interaction between pairs 
of parallel bonds, is sufficient to transform the melting transition from second- 
order, as seen in the case a = b = (original Flory model), to first-order. If 
b is negative, the melting and mode-coupling temperatures decrease. If b is 
positive, instead, these temperatures increase and the melting transition turns 
into a second-order transition. This is true for any positive value of b in the case 
a = 0, while it is true for large enough negative values of b when a is positive. 

It is interesting to observe that the parallel bond density in ML and in CR 
is always unity while its value for the EL depends on the temperature. 

4.1.3 a>l 

We consider now the case a > 1. This case corresponds to a ground state that 
is not crystalline, as shown in Fig. 4(b). When the four-site interaction is 
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Figure 11: a = 0.5,5 = 0. Specific heat in the 2-cyclc FP scheme for the ML 
(continuous hne), EL/SCL (dashed hne) and CR (dash-dot line). We also show 
Tcre(^), Tm(o) and Tmc(»). 



stronger than the three-site interaction, the polymer assumes a configuration 
that is such that the number of parallel bonds is minimized. In this case there 
is a very high number of gauche conformations at low temperatures and, even 
though the polymer assumes an ordered configuration on the lattice, it is not 
a crystalline configuration according to our definition. Therefore, we do not 
consider this case any further. 

4.1.4 a<0 

It is also possible to analyze the case in which a < 0. When a is negative 

(corresponding to a negative four-site interaction energy Sp) the temperature 
Tmc at which the EL appears moves to higher values. Since the temperature 
Tm(= Tore) at which CR appears is unaff'ected by the choice of the value of 
a, a shift of the origin of the EL phase to higher temperatures makes ML the 
equilibrium phase for temperatures between Tm and Tmc > Tm; it is no longer 
a metastable phase in this range. We identify the equilibrium portion of the ML 
phase as a new equilibrium phase, and denote it by ELml! the subscript is a 
reminder that the phase is associated with the ML phase. The new phase is 
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Figure 12: (a) Gauche bond and (b) parallel bond density in the 2-cycle FP 
scheme for the ML (continuous line), EL/SCL (dashed line) and CR (dash-dot 
line). We also show Tcre(^), T'm(o) and Tmc(»)- 



again a liquid phase. Hence, the transition at Tmc is a liquid -liquid transition 

and is continuous. Similarly, the transition at Tm between CR and ELml is 
also a continuous transition. Since ML exists below Tm, we can formally treat 
the ML phase below Tm as obtained by continuation of ELml below Tm, but 
this continuation should not be taken as a supercooled liquid below Tm, as there 
will be no energy barrier due to the continuous melting transition. Furthermore, 
since the liquid-liquid transition at Tmc occurs at a temperature higher than the 
melting tcmpcrat\.irc, this case has no relevance for st\.idying the glass transition. 
Hence, we do not pursue it further. The same behavior is observed, as explained 
above, when a = and 6 > 0. 

4.1.5 Choice of a 

Because of the above considerations, we consider a in the range <a <1. In 
this range, the four-site interaction is repulsive (the system spends some en- 
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Figure 13: Effect of b on the phase diagram of the system (a = 0.5). The free 
energy of ML and CR does not depend on b. Three EL/SCL curves are shown 
corresponding to b = —0.4 (long dash), b = (medium dash) and b = 0.4 (short 

dash). 



ergy to ahgn two segments parallel to each other) and the ground state is the 
crystalline one [see Fig. 4(a)]. In this case, the model exhibits a first order 
melting transition at a temperature Tm = 7m (a) between EL, that is stable at 
temperatures higher than Tm, and CR, that is stable for temperatures lower 
than Tm- It can be observed that the discontinuity in the specific heat at Tmc 
is a function of the parameters a and b. In particular, if we fix b, as a increases 
the discontinuity gets smaller and smaller as long as a <0.8, and then starts 
growing again while the transition temperature keeps moving to lower values; 
the results are not shown. 

The crystalline phase is an ordered one but, unlike the ground state predicted 
by the original Flory model [1,2] (Fig. 1), it has non zero entropy. It also 
satisfies the Gujrati-Goldstein bounds. The I and O states disappear in the 
crystalline phase but this phase has non zero entropy because of the many 
possible configurations that the polymer chain can assume corresponding to 
different sequences of L and R states. The entropy of the crystalline phase goes 
to zero only when the temperature goes to zero, which is consistent with the 
third law of thermodynamics. 

If the cooling process is such that the system can avoid crystallization at Tm, 
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the equilibrium liquid EL can be supercooled to give rise to SCL that transforms 
into the metastable liquid ML through a liquid-liquid second-order transition at 
Tmc < (no latent heat is associated with the transition). The metastable 
liquid and the equilibrium liquid phases arc somehow similar. The metastable 
liquid consists of a random sequence of R and L states, while the equilibrium 
liquid consists of a random sequence of R, L, I and O states. The presence of 
O and I states makes the total energy and the entropy of the equilibrium liquid 
larger than those of the metastable liquid, see Figs. 10 and 11. 

4.2 Relation with the Mode-coupUng Transition 

We tentatively identify the critical temperature Tmc of the liquid-liquid transi- 
tion as the mode-coupling transition temperature because the transition exhibits 
some of the features predicted by the original mode coupling theory at the criti- 
cal mode-coupling transition temperature. It should be remarked that our equi- 
librium investigation cannot provide any direct information about the dynamics 
at Tmc except by inference. Hence, the connection we allude to above should 
only be taken as tentative, in view of the fact that the mode-coupling transition, 
applied successfully to simple fluids, is considered to be a dynamic transition. 
We can only add that the mode-coupling theory is not well-understood for long 
polymers, and it is not clear what its predictions might be for infinite polymers 
that we are investigating here. 

According to this theory, the dynamics slows down according to Eq. (3) 
near Tmc- The local molecular structure freezes and only long-time cooperative 
jumps are allowed below this temperature. Thus, the dynamic transition is 
between two disordered states, very much like the thermodynamic liquid-liquid 
transition we observe in our calculation. Let us consider the behavior of the 
correlation length ^scl of the system near the critical temperature Tmc- As 
Tmc is approached from above (T T^q), the correlation length Cscl of the 
supercooled liquid must diverge to infinity because the transition is continuous. 
It is very easy to observe from the results that there is a discontinuity in the 
specific heat of the system at the transition from the supercooled liquid to 
the metastable liquid. Indeed, the SCL terminates at Tmc as T ^ ^mc- '^^'^ 
disappearance of SCL is what gives rise to this divergence, which will contribute 
to the critical slowing down of the system. Such a critical slowing down is 
exhibited at the mode-coupling transition [39]; see Eq. (3). 

On the other hand, ML exists at all temperatures. Thus, there would be no 
divergence at Tmc in the correlation length ^ml associated with ML . Indeed, 
its specific heat remains continuous. This will imply that the dynamics of the 
system should not undergo any significant change at the critical temperature 
Tmc if we approach it while heating up ML in such a way that ML is not allowed 
to turn into SCL. In simulation, one can investigate ML by suppressing fast 
relaxations that are supposed to freeze at the mode-coupling temperature. Such 
an attempt has already been made [25] where one sees no anomalous behavior 
at the mode-coupling temperature. Parisi and coworkers [25] while analyzing a 
Lennard-Jones system have observed this kind of dynamics. In their approach. 
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the fast dynamics of the system (the one pertinent to the supereooled Uquid 
in our model) is suppressed and only the slow dynamics of the system is taken 
into account. The slow dynamics is described as a relaxation process taking 
place in a connected network of potential energy minima. Indeed, the authors 
only observes an Arrhenius behavior in the relaxation time, as opposed to the 
Vogel-Tammann-Fulcher behavior associated with the mode-coupling transition. 
Even though the system studied by Parisi and coworkers is very different from 
the polymer system studied here, it is important to note that all the numerical 
results obtained in the case of the Lennard- Jones fluid are in agreement with the 
experimental findings in non-network forming glasses and especially in glasses 
that are fragile according to Angell's definition [40] . 

Recent experimental results obtained by Sokolov and coworkers [41] studying 
polyisobutylene and polystyrene show the presence of a critical behavior only 
above Tmc along with the failure of the predictions of the mode coupling theory 
below the critical temperature. These authors show many similarities between 
the results obtained on these polymeric glasses and the ideas of the liquid-liquid 
transition in polymeric liquids formulated by Boyer and coworkers [42,43]. It 
is worth noting that this liquid-liquid transition would manifest itself through 
a discontinuity in the first derivative of the specific heat (and not the specific 
heat itself as in the present case) at the transition temperature. This makes 
Boyer's result very different from our result. The idea introduced by Boyer has 
been strongly criticized and is still the subject of discussion [44,45]. 

The second similarity has to do with our choice of the parameter a, so that 
Tmc lies above the Kauzmann temperature Tk- This is also what is expected in 
the mode-coupling theory in which the transition occurs above the conventional 
glass transition temperature. 

The third similarity appears when we allow free volume in our model in Eq. 
9(b), as has been done recently [29]. It is found that the free volume in the 
model for the case of infinite polymers vanishes identically at TmC) and remains 
zero below it. Consequently, one expects the viscosity to diverge at Tmc- 

While the mode coupling theory describes the transition at Tmc as dynamic 
in nature, our results show that the transition at Tmc is thermodynamic in 
nature. The sharp transition is due to the polymer being infinitely long, and 
disappears as soon as polymers become finite in size [29]. However, for polymers 
of reasonable sizes, there will continue to be a narrow crossover region between 
two phases (ML and SCL). 

4.3 Ground State and Kauzmann Temperature 

Below Tmc, the only two phases that are present are the metastable liquid and 
the crystal. Above Tmc, the supercooled liquid, which is the continuation of 
the equilibrium liquid above Tm, is more stable than the metastable liquid and 
coexists with the crystal. It is worth noting that the modified free energies F of 
both the metastable liquid and the supercooled liquid cross over and becomes 
positive at some finite and non-zero temperature. Let us focus on the metastable 
liquid as its behavior is easy to describe since its modified free energy F re- 
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mains independent of a and b. Wc first observe that the 2-cyclc FP solution 
contains within its possible solutions the 1-cycle solution. We also find that the 
free energy of all possible 2-cycle solutions (including the 1-cycle solutions) at 
absolute zero arc the same: F = 0. 

Because of the exact nature of our calculation, this equality of the ML and 
CR free energies at absolute zero is not brought about due to any approxi- 
mation. Because of this equality at absolute zero, we will now consider the 
modified free energy F in the following. The CR free energy remains negative 
at all temperatures and approaches zero at absolute zero. JThus, CR has non- 
negative entropy. On the other hand, the ML free energy F, which is negative 
at higher temperatures, becomes positive at some intermediate and non-zero 
temperature Tmo and keeps increasing, as the temperature is lowered, until the 
Kauzmann temperature Tk (~ 0.48 in our model) is reached. Below Tk, the 
ML free energy F must jiecessarily decrease since it must vanish at absolute 
zero. The maximum of F corresponds to the vanishing of the entropy of the 
ML phase, below which the entropy must become negative [46]. Thus, the ex- 
istence of the^Kauzmann temperature is a consequence of the fact that the ML 
free energy F , once it crosses the zero at Tmo, must necessarily decrease at 
some lower temperature so as to return to zero at absolute zero. The existence 
of the maximum in Fml as a function of the temperature is the root of the 
rapid entropy drop noted by Kauzmann [5] . The maximum at a positive Tk is 
forced by thermodynamics since the larger specific heat of ML makes Fml cross 
over to positive values at Tmo- If we had carried out our calculation in some 
approximation, as is the case with the calculation of Gibbs and DiMarzio [8], 
we certainly could not draw this remarkable conclusion. 

4.4 S'cR>5'scL and Entropy Crisis 

The crystalline phase has an entropy that is never negative. Hence, its entropy 
is larger than the entropy of the metastable liquid in the temperature interval 
Tk < T < Toq, where Teq is the temperature at which the entropy of the two 
phases is the same (see Fig. 10(b)). This result contrasts the common belief [17] 
that the entropy of a crystalline phase must always be lower than the entropy of 
the corresponding liquid phase, even if the latter is an equilibrium phase. Our 
results clearly show that there is no thermodynamic requirement for this belief 
to be true. Indeed, real systems like He conform with this observation. 

In order to sustain the common belief that the entropy of the liquid must 
always be larger than that of the crystal, it was conjectured by Kauzmann 
that the system must avoid the (Kauzmann) catastrophe caused as soon as the 
requirement Sqr < S'scL is violated. The system is supposed to avoid the 
catastrophe by undergoing either a spontaneous crystallization, as proposed by 
Kauzmann in his original paper, or an ideal glass transition [5,8,15,47]. 

The most important result of the present research regarding the Kauzmann 
catastrophe is that since our calculations show that it is possible to observe 
'S'cR > 'S'scL in an explicit model calculation, the existence of the catastrophe 
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must be reinterpreted in terms of the entropy crisis corresponding to having no 
reahzable state (negative entropy) [46]: the catastrophe happens at a tempera- 
ture Tk where the entropy of the metastable Hquid goes to zero and not at the 
temperature Tcq at which the entropy of the hquid becomes equal to the entropy 
of the crystaL The two temperatures coincide in the original Flory model since 
the entropy of the crystal is identically equal to zero below Tm, see Fig. 1, while 
in our exact calculation, the two temperatures are different because the entropy 
of the crystal is zero only when the temperature goes to zero. 

4.5 Ideal Glass Transition 

The ML free energy has no physical relevance below Tk since it corresponds 
to negative values of the entropy. Decreasing the temperature as T ^ the 
entropy falls very rapidly to zero (as seen, for example, in Fig. 10 (b)). As the 
temperature is decreased, the energy of the metastable liquid cannot increase 
because this would correspond to a negative specific heat. At the same time, 
this energy cannot decrease since there arc no states with non-negative entropy 
(except the crystal phase, but wc do not consider this possibility here) available 
to the system. The only conclusion that can be drawn from these observations 
is that for T < Tk, the metastable liquid remains frozen in the state in which it 
finds itself at Tk- This describes the ideal glass transition. In our analysis, ML 
does not undergo any changes in its state at the glass transition, so that this 
transition cannot be a first-order transition as recently proposed by Parisi and 
coworkers [26(b)]. 

It is interesting to notice that the energy of the metastable liquid increases 
monotonically with T and is very similar to the excitation profiles for other 
systems like the mixed Lennard-Jones system [48] and the two states model 

[49]. .... 

It is also interesting to notice that our model predicts, at the Kauzmann 
temperature, an upward discontinuity for the specific heat. The specific heat of 
the metastable liquid is decreasing with the temperature for T > Tk- This kind 
of behavior is in disagreement with many experimental results but has been ob- 
served in computer simulations by Parisi and coworkers [26] and experimentally 
(at least in some temperature interval immediately above the glass transition) 
in glasses formed by low molecular weight materials of very different nature like 
1-butene [45] and the metallic system Au-Si [40] . 

4.6 Flory Model 

The value a = 1 corresponds to a borderline case. As we have shown in Fig. 4, 
if a is larger than 1 then the ground state is not crystalline anymore and there 
arc no parallel bonds at T = 0. Correspondingly, we observe from our results 
that as a — *■ 1, Tmc moves to lower and lower values and eventually goes to zero 
when a goes to one. We also find that, if we keep b fixed, then as the value of 
a increases, the values of Tm and Tmc decrease and, as said above, for a > 0.8, 
the supercooled liquid shows its own Kauzmann temperature corresponding to 
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Figure 14: Dependence of the melting temperature on o. The first-order tran- 
sition fine (continuous Hne ) and the second-order transition Une (dashed Une) 
are shown. 



the maximum in its free energy. If a becomes negative, the critical temperature 
moves to higher temperatures. Unlike the case of o > 0, there is a temperature 
interval (Tm < T < Tmc) where the mctastable liquid becomes the stable phase 
ELml- The melting transition that is observed is a continuous transition. This 
aspect is in disagreement with experimental observations. 

If a = 6 = 0, as explained before, the model reduces to the Flory model of 
polymer melting. In this case Tm = Tmc = Tore and the melting transition 
from the equilibrium liquid to the crystalline phase becomes continuous in con- 
trast with the original calculations of Flory. It is important to notice that in this 
case there is no possibility to obtain a supercooled liquid since the equilibrium 
liquid phase disappears continuously into the crystal at Tm- 

The solution foimd by Flory for the original model corresponds to a tricrit- 
ical solution: if we consider the melting transition present in the system, this 
transition is first order for a > and continuous for a < 0. In the latter case, 
wc have Tm = Tcre, as shown in Fig. 14. 

It seems reasonable to assume that in order to be able to describe the physics 
of real systems, the value of the parameter a must be chosen in the range between 
and 0.7 while the value of b should be in the range —0.5 — 0.5. 
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4.7 Compcirison with A Real System 

By a proper choice of various parameters in our model, we can fit the predictions 
of our theory with experiments. We discuss one such example below. Setting 
a = 0.5 and 6 = 0, for example, our model predicts Tyi/T^^ ~ 1.20/0.48. We can 
consider polyethylene (PE) and try to describe its thermodynamic properties 
using our model. The experimentally measured melting temperature of PE is 
Tm(PE) ~ 400K. Then the model predicts Tk(PE) ~ 160K. This temperature 
is about 40K below the experimentally determined glass transition. Since we 
expect the experimental glass transition to occur above the ideal glass transition 
because of experimental constraints, we can conclude that the prediction of our 
model is, at least, reasonable. 

5 CONCLUSIONS 

We have considered an extension of the Flory model of melting by introducing 
two additional interactions characterized by parameters a and h. One interaction 
is between a pair of parallel bonds and the second one is due to a hairpin turn. 
The model is defined on a checkerboard version of the square lattice, and has 
been solved exactly on a Husimi cactus, which is a recursive lattice [36]. It 
should be recalled [36] that calculations done on a recursive lattice have been 
shown to be more reliable than conventional mean- field calculations. The choice 
of the Husimi cactus is also important for the inclusion of the Gujrati-Goldstein 
excitations that are responsible for destroying the complete order in the crystal 
phase CR [10-12] in the Flory model. The method of calculation is to look for 
the fix-point (FP) solution of the recursion relations. We need to consider 1- 
cycle and 2-cycle FP solutions to describe the disordered phase and the crystal 
phase, respectively. This has required us to provide in this work the extension 
of the Gujrati trick [36] to calculate the free energy of the 1-cycle FP solution 
to the 2-cycle solution. The exact nature of the solution allows us to draw 
some important conclusions, which would not have been drawn with the same 
force, had we obtained the solution under ambiguous and/or uncontrollable 
approximations. 

We have identified an equilibrium liquid phase EL at high temperatures 
T > Tm- There is another liquid phase ML which exists at all temperatures, 
but which never becomes an equilibrium phase for appropriate choices of the 
parameters a and b. Below the melting temperature Tm, the crystal phase CR 
becomes the equilibrium phase. The melting transition is usually a first order 
transition with a latent heat, below which we can continue EL to give rise to 
the supercooled liquid phase SCL. This phase terminates at a lower temperature 
Tmc, where it meets ML continuously with no latent heat. The transition at 
Tmc is a continuous liquid-liquid transition. 

Both SCL and ML represent mctastable phases in the system. For a metastable 
state to exist in Nature, its entropy must be non-negative. A negative entropy 
{S < 0) in metastable states (SCL and/or ML) implies that such states cannot 
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exist. Wc have argued that the ML free energy must have a maximum at Tk 
as a function of the temperature T; see Fig. 8. Thus, ML has non-negative 
entropy above Tk, but gives rise to negative entropy below Tk- We have called 
the appearance of 5* < the entropy crisis , which we have used instead of the 
Kauzmann paradox (5ml < Scr) as the driving force behind the ideal glass 
transition at Tk- The ideal glass transition occurs in the metastable liquid ML, 
and not in SCL, which is contrary to the conventional wisdom. The portion of 
ML below Tk must be replaced by the ideal glass (see dotted thin horizontal 
portion below Tk in Fig. 8), which is completely inactive in that its entropy and 
specific heat are both zero. The rapid drop in the entropy near the Kauzmann 
temperature Tk is a direct consequence of the existence of the maximum in the 
ML free energy. The liquid-liquid transition at Tmc between the two metastable 
phases ML and SCL has been shown to share many similarities with the critical 
mode-coupling transition, even though the latter is known to be driven by dy- 
namics in simple fluids. It should be remarked that nothing is known about this 
dynamic transition in the semiflexible Hamilton walk limit of infinite polymers. 

The Flory model is shown to give rise to a continuous melting. Indeed, 
the melting point in the Flory model turns out to be a tricritical point in our 
calculation. It would be interesting to see if this conclusion can be sustained in 
other computational scheme. 

The natural extension of this model in Eq. (9b) involves the analysis of the 
effects of compressibility (taking into account voids as another species on the 
lattice) and of finite chain size, both in the polydisperse and in the monodisperse 
case. This is reported elsewhere [29,30]. 

Wc finally discuss an interesting aspect of the 2-cycle FP solution for 
EL/SCL observed by Semerianov [50]. The values of ia, ih, Oa, and Ob de- 
pend on the choice of the initial guesses used in the FP solution of the recursion 
relations. Thus, there are many different 2-cycle solutions for EL/SCL that dif- 
fer in the values of i^,, i^,, Oa, and Ob.What we find is that the product io remains 
the same on both sublattices for all initial guesses, and that they all give the 
same free energy and densities. We hope that the observation will provide some 
useful information about the free energy landscape picture, but this requires 
further investigation. We hope to report on this in near future. 
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APPENDIX 



We index the cactus levels in a slightly different manner here for simplicity. 
If the base site of a square is indexed m, then the intermediate and the peak 
sites are all indexed m+1. The square is still called the mth square. The origin 
is indexed, as before, m = 0. Wc evaluate Af(m), the number of sites at the mth 
generation of a tree rooted at the origin (m = 0). The rooted tree is only half 
of a complete tree. It is evidently given by: 

so that the total number of sites belonging to the first m generations (and 
excluding the origin) of a rooted tree is: 

fe=i ^ 

The total number of sites belonging to the first m generations of the complete 

tree is then equal to twice the total number of sites belonging to the first m 
generations of a rooted tree plus one since wc have to consider the origin too: 

Nm = 2A^(™^*"'' + 1 = 3™+! - 2. 

Let us consider now the number of the mth generation squares (m > 0) 
of the rooted tree. Clearly we have: 

C nm 

'^(m) — "J • 

The total number of squares belonging to the first m generations of a rooted 
tree (so that the maximum generation of the squares is m — 1) is: 

om 1 

arooted \ ^ a 

•^M -Z^^ik)- — 2 — ■ 

fc=0 

The total number of squares belonging to the first m generations of the complete 
tree is then just twice the total number of squares belonging to the first m 
generations of a rooted tree: 

Q o QTOoted om -i 

In order to make the cactus homogeneous, jwe must consider it to be part of 
a larger cactus. Consequently, both Sm and must be modified in order to 
take into account the presence of dangling bonds at surface sites: we consider 
to add half a square to each surface site. Recalling that each square contains 
four half-sites (a site is shared by two squares), we conclude that each square 
contributes two sites to the number of sites. Hence, each half-square contributes 
one site to the site count. Thus, we modify Nm by adding 1 site for each surface 
square. This gives for the total number of sites 

A^m = A^m + 3™ = 4 X 3"^ - 2. 
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Modifying Sm~i by the half-squares at the surface sites gives for the total 
number of squares: 

Sm = Sm-l + ^ X 2 X 7V(„) = 2 X 3™ - 1. 
Now, if we consider the thermodynamic limit in which m — > oo we have: 




which is consistent with our earlier homogeneous hypothesis Ns = N relating 
the total number of squares TVs f^ud the total number of sites N. 

Let us finally consider the total number of bonds in the first m gene£a- 
tions (m > 0) of the tree. Each square contains four bonds, and there are Sm 
squares in this tree. Thus _ 

Bm = 4(3" - 1). 

The modification of the lattice introduced above implies that we must now add 
half-square at each of the surface sites of the mth generations (m > 0) tree. 
Each half-square contributes 2 bonds. Hence 

B„ = 5„ -h 4 X 3™ = 8 X 3" - 2. 

If we consider the thermodynamic limit in which m — > oo we have: 




as expected. 
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